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BLIND CARRIER OFFSET DETECTION FOR QUADRATURE MODULATED 
DIGITAL COMMUNICATION SYSTEMS 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention relates generally to digital communication systems using quadrature 
modulation techniques. More specifically, the invention relates to a system and method for blind 
detection of carrier frequency offsets in such systems. 

Description of the Prior Art 

A digital communication system typically transmits information or data using a continuous 
frequency carrier with modulation techniques that vary its amplitude, frequency or phase. After 
modulation, the signal is transmitted over a communication medium. The communication media 
may be guided or unguided, comprising copper, optical fiber or air and is commonly referred to 
as the communication channel. 

The information to be transmitted is input in the form of abit stream which is mapped onto 
a predetermined constellation that defines the modulation scheme. The mapping of each bit as 
symbols is referred to as modulation. 

Each symbol transmitted in a symbol duration represents a unique waveform. The symbol 
rate or simply the rate of the system is the rate at which symbols are transmitted over the 
communication channel. A prior art digital communication system is shown in FIG. 1 . While the 
communication system shown in FIG. 1 shows a single, communication link, those skilled in this 
art recognize that a plurality of multiple access protocols exist. Protocols such as frequency 
division multiple access (FDMA), time division multiple access (TDMA), carrier sense multiple 
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access (CSMA), code division multiple access (CDMA) and many others allow access to the same 
communication channel for more than one user. These techniques can be mixed together creating 
hybrid varieties of multiple access schemes such as time division duplex (TDD). The type of 
access protocol chosen is independent of the modulation type. 

One family of modulation techniques is known as quadrature modulation and is based on 
two distinct waveforms that are orthogonal to each other. If two waveforms are transmitted 
simultaneously and do not interfere with each other, they are orthogonal. Two waveforms 
generally used for quadrature modulation are sine and cosine waveforms at the same frequency. 
The waveforms are defined as 

s { (t) = y4cos(2 7rfJ) Equation 1 

and 

s 2 (t) = A sin(2 Ttfj) Equation 2 

where f c is the carrier frequency of the modulated signal and A is the amplitude applied to both 
signals. The value of A is irrelevant to the operation of the system and is omitted in the discussion 
that follows. Each symbol in the modulation alphabet are linear combinations generated from the 
two basic waveforms and are of the form a x cos(2 7f c t) + a 2 sin(2 7fj) where a, and a 2 are real 
numbers. The symbols can be represented as complex numbers, a } + ja 2 , where j is defined as j 
= 7-1. 

The waveforms of Equations 1 and 2 are the most common since all passband transmission 
systems, whether analog or digital, modulate the two waveforms with the original baseband data 
signal. Quadrature modulation schemes comprise various pulse amplitude modulation (PAM) 
schemes (where only one of the two basic waveforms is used), quadrature amplitude modulation 
(QAM) schemes, phase shift keying (PSK) modulation schemes, and others. 
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A prior art quadrature modulator is shown in FIG. 2. The modulator maps the input data 
as a pair of numbers {a l9 a 2 } which belong to a set defined by the modulation alphabet. a x 
represents the magnitude (scaling) of the first waveform and a 2 represents the magnitude (scaling) 
of the second waveform. Each magnitude is modulated (i.e. multiplied) by the orthogonal 
waveforms. Each individual modulator accepts two signal inputs and forms an output signal at the 
carrier frequency. 

A prior art quadrature demodulator is shown in FIG. 3. The demodulator generates sine 
and cosine waves at a carrier frequency f c for demodulation. Ignoring channel effects, the received 
signal can be represented as 

r(t) = a l (t)cos(2nfj + ^ Q ) + a 2 (t)sin(2nf c t + ^ Q ) Equation3 

where a x {i) represents the plurality of amplitudes modulated on waveform s\(t) as defined by 
Equation 1 and a 2 (t) represents the plurality of amplitudes modulated on waveform s 2 (t) as defined 
by Equation 2. $ is an arbitrary phase offset which occurs during transmission. 
The cosine and sine demodulator signal components are defined as: 

r c (0 = r(0*cos(2^0 

= i a x COS((/ c - f w )t + <f> 0 ) + ^a 2 Sin((/ C - f LO )t + 0 o ) Equation 4. 

+ cos((/ c +f w )t + </> 0 ) + ^a 2 sm((f c -f w )t + 0 o ) 



and 
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r,(0 = KO*sin(2^0 

= ^a 2 cos((/ c -f w )t + ^)-^ a i sm((/ c -/ioX + ^ 0 ) Equation5 

-^fl 2 cos((/ c +/ iO > + ^ 0 ) + ^i sin((/ c -f LO )t + fa) 



The carrier frequency components,^ + f W9 are suppressed by the lowpass filters. The 
signals after filtering are: 



y c it) = \a, cos((/ c -f w )t + <fi 0 )+±a 2S m((f c Equation6 



and 



^(0 = ^ 2 cos((/ c -f w ) t+ t 0 )-L ai sin((/ c + Bqwltai7 



If the local oscillator frequency in Equations 6 and 7 is equal to the carrier frequency,^ 
=f a and the phase offset is equal to zero, <j> 0 = 0, the right hand sides of Equations 6 and 7 become 
Vi a x (t) and V2 a 2 (t) respectively. Therefore, to effect precise demodulation, the local oscillator 
must have the same frequency and phase as that of the carrier waveform. However, signal 
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perturbations occurring during transmission as well as frequency alignment errors between the 
local oscillators of the transmitter and receiver manifest a difference between the carrier and local 
oscillator frequencies which is known as carrier offset. A phase difference between the carrier 
and local oscillator frequency is created as well. However, if the difference in frequencies is 
corrected, the difference in phase is simple to remedy. Phase correction is beyond the scope of 
the present disclosure. 

Carrier frequency offset is defined as: 

A/ = f c - f LQ . Equation 8 

To synchronize either parameter, the frequency and phase offsets need to be estimated. In 
prior art receivers, frequency offset estimation is performed after a significant amount of data 
processing. Without correcting offset first, the quality of downstream signal processing suffers. 

What is needed is a system and method of detecting and estimating carrier frequency offset 
before any data signal processing is performed. 

SUMMARY OF THE INVENTION 
The present invention provides a system and method for estimating the frequency offset 
experienced between carrier and local oscillator frequencies in communication systems using 
quadrature modulation. The invention exploits the geometry of the quadrature modulation 
constellation and estimates actual offset within a predefined carrier offset value without requiring 
data estimation. 

Accordingly, it is an object of the invention to provide a less-complex system and method 
for blindly estimating carrier frequency offset. 

It is another object of the present invention to blindly estimate carrier offset in a 
communication system using quadrature modulation regardless of the access protocol. 
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Other objects and advantages of the system and method will become apparent to those 
skilled in the art after reading a detailed description of the preferred embodiment. 



BRIEF DESCRIPTION OF THE FIGURES 
FIG. 1 is a simplified system diagram of a prior art digital communication system. 
FIG. 2 is a system diagram of the prior art quadrature transmitter shown in FIG. 1. 
FIG. 3 is a system diagram of the prior art quadrature receiver shown in FIG. 1. 
FIG. 4 is a system diagram of the blind carrier offset estimator of the present invention. 
FIG. 5 is a detailed system diagram of a blind digital carrier offset estimator of the present 
invention. 

DETAILED DESCRIPTION OF THE INVENTION 
The embodiments will be described with reference to the drawing figures where like 
numerals represent like elements throughout. 

Shown in FIG. 4 is an analog or digital blind carrier detector 33 of the present invention. 
A quadrature modulated signal r(t) is received from a communication channel (not shown) and is 
input 19 to a receiver 17. One skilled in this art recognizes that additional conversion means may 
exist before the detector input 19 to convert the energy used in the transmission media to 
compatible signals and is beyond the scope of this disclosure. The received signal r(i) is coupled 
to a cosine mixer 21 c and a sine mixer 21 s . Each mixer 21 c , 21 s has a first input 25 c5 25, for 
coupling with the received signal r(t) and a second input 27 w 27 5 for coupling with the output of 
a local oscillator LO. The local oscillator LO is programmed to generate cosine and sine waves 
at the carrier frequency./; (Equations 4 and 5) of the received signal r(i). 
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The carrier-frequency demodulated outputs r c (t), r 5 (0 &om each mixer 21 c , 21 s are input 
to respective lo wpass filters 29^ 29 8 which suppress high-frequency noise components impressed 
upon the received signal r(t) during transmission through the transmission media and mixer sum 
frequencies, f c + f W9 (Equations 6 and 7). As in prior art demodulators, the response 
characteristics of the lowpass filters 29<, 29 a may be a bandwidth as narrow as AfMAx - the 
maximum allowable carrier offset. The output y c ((), y s (t) from each lowpass filter 29 c , 29, is 
coupled to inputs 31c 31 s of a carrier offset estimator 33. 

The carrier offset estimator 33 produces an estimate of the carrier offset 35 before data 
signal processing commences using a complex power processor 37 in conjunction with a complex 
Fourier transform processor 39. The filtered, carrier frequency demodulated cosine and sine 
components of the quadrature signal y c (t) an.&y s (t) are coupled to the complex power processor 
37 which performs an intermediate power calculation of each quadrature component in the form 
of where the powers j> comprise integer multiples of four; i.e. y - 4, 8, 12, 16 .... In the 
preferred embodiment, the powers is 4. 

The complex power processor 37 may be implemented to raise the input complex signal 
to a power which is any positive integer multiple of four. Caixier offset detection systems which 
use a complex power processor with a power of two or its positive integer multiples are known 
in the art. However, these prior art systems do not work in quadrature-modulated digital 
communication systems. To properly detect a carrier offset in a quadrature-modulated digital 
communication system demodulator, a complex power of four or its integer multiples are 
necessary. 

The complex power processor 37 combines the lowpass filter outputs^) and y s (f) into 
a single complex value signal y(i) defined as: 
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y{t) =y c (t) +M0 Equation 9 

where j is defined as 7 = \M. The complex power processor 37 generates two power output 
signals 

q c (t) = Re{(y(f)y} Equation 10 

and 

q s (t) = Im {(KO) 4 } Equation 1 1 

where Re{x) denotes the real part of a complex number x 9 and Im{x} denotes the imaginary part 
of the complex number x. The complex power processor 37 removes the modulation component 
from each received symbol leaving the carrier frequency. The real q c {i) and imaginary q s {t) signal 
components are output and coupled to the complex Fourier transform processor 39. 

The complex Fourier transform processor 39 treats the real q c {() and imaginary q s (t) signal 
components as a single complex input signal q(i) = q c (i) + jq s (f). The processor observes q(t) for 
a finite period of time TV and computes a complex Fourier transform of the observed signal q(t) 
over this period of time. 

The Fourier processor 39 performs a Fourier transform of the power processed signals 
from the observed period J^and outputs a frequency at which the amplitude of the transform was 
measured to be maximal during that time period T w . The output 35 represents an accurate 
estimate of bf and is signed since the transform input signal is complex. The sign identifies 
whether the local oscillator LO frequency is less than or greater than the carrier frequency. 

A detailed, low-complexity digital implementation of the present invention 53 is shown 
in FIG. 5. Lowpass filter 29 w 29 8 output signals y c {i) and y s (i) are sampled at a sampling rate/ 
to produce discrete-time signals y c [n) andy s [n]. To ensure that all possible carrier frequency 
offsets up to A/max are detected, 24^^ < / 8 must be satisfied. The passband of the low pass 

-8- 
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filters 29 c , 29 s must be wider than /Vmax t0 avoid suppressing the signal which contains the carrier 
offset information. 

The sampled signals y c [n) and y s [n] are input 51 c , 51 s to a complex power processor 57 
and combined as a single complex signal, y[ri\, where y[n] -y c [n] + jy s [n\. The power processor 
57 produces a complex output defined by q[n] = (y[«]) 4 . The output q[n] is coupled to a buffer 
59 for accumulating N outputs from the complex power processor 57. 

The accumulated block of complex numbers Nis coupled to a digital Fourier transform 
(DFT) processor 61 which performs a transform from the time domain to the frequency domain for 
the ^complex numbers. The DFT processor 61 outputs Ncomplex numbers corresponding with 
the input N. Each number is associated with a particular frequency ranging from -ffl to {+fJ2 - 
fJN). Each frequency is^i/ATaway from a neighboring frequency. The frequency domain values 
outputbytheDFT61 are assembled and compared with one another. The value having the largest 
magnitude represents the best estimate of the carrier frequency offset Af. 

The embodiment described in FIG. 5 is capable of estimating all carrier frequency offsets 
smaller than^/2. This follows from the restriction 24^^ < f s imposed above. The carrier offset 
is resolved to within a frequency uncertainty of ± f/lN since the frequencies at the output of 
the DFT 61 are quantized to a grid with a spacing of fs/N. Since the frequencies are fsIN away 
from each other, the invention 53 renders precision within ±^ of the selected value. Therefore, 
the number of samples N accumulated for the Fourier processor 61 to transform determines the 
resolution of the carrier offset estimate Af. An efficient implementation of the DFT 61 used in the 
present invention 53 can be achieved using the fast Fourier transforms (FFT) family of algorithms. 

The present invention 33, 53 may be physically realized as digital hardware or as 
software. The lowpass filters shown in FIG. 5 may be realized in digital hardware or software 
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operating at a sampling rate faster than f s . In some communication systems, for example those 

employing CDMA protocols, the lowpass filters and downsamplers f s may be replaced with 

accumulators and integrate-and-dump processes. 

While the present invention has been described in terms of the preferred embodiments, 

other variations which are within the scope of the invention as outlined in the claims below will 

be apparent to those skilled in the art. 

* * * 
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CLAIMS 

What is claimed is: 

1 . A detector (33) for use in a receiver (17) that estimates a difference (A/) between a 
received quadrature demodulated communication signal (r(t)) carrier frequency (f c ) and a local 
oscillator (LO) frequency used for demodulation, the detector (33) characterized by: 

a complex power processor (3 7) for receiving the quadrature demodulated communication 
signal (r(f)) cosine (y c (t)) and sine (y s (t)) components producing real (Re{(y(t)y}) and imaginary 
(Im{(y(t)Y}) power signal components, said power signal components coupled to a complex 
Fourier transform (39); 

said complex Fourier transform (39) performing a complex Fourier transform on said 
power signal components for a predefined period of time outputting a frequency transform 
containing a plurality of frequencies observed during said predefined period of time; and 

said plurality of frequencies having a maximum frequency corresponding to the difference 
(Af) between the carrier (Q and local oscillator (LO) frequencies. 

2. The detector (33) according to claim 1, wherein said complex Fourier transform (39) 
maximum frequency is further characterized with a sign (±) indicating whether the difference 
frequency (A/) is greater than or less than the local oscillator (LO) frequency. 

3. The detector (33) according to claim 2, wherein said complex power processor (37) 
is further characterized by obtaining a power f) that is an integer multiple of four. 
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4. A method (33) for estimating a difference (A/) between a received (17) quadrature 
demodulated communication signal (r(0) carrier frequency (f c ) and a local oscillator (LO) 
frequency used for demodulation, the method (33) characterized by the steps of: 

a) performing a complex power calculation on the received quadrature demodulated 
communication signal (r(f)) cosine (y c (t)) and sine (y/f)) components producing real (Re {(y(t)Y}) 
and imaginary (Im{(y(t)y}) power signal components; 

b) transforming said real (Re{(y(f)Y}) and imaginary (Im{(y(t)Y}) power signal 
components into Fourier frequency domain values representing a plurality of frequencies; 

c) selecting a maximum frequency from said plurality of frequencies; and 

d) outputting said selected frequency as the difference frequency (A/). 

5. The method (33) according to claim 4 wherein step d) is further characterized by a step 
of assigning a sign (±) to said maximum frequency indicating whether the offset frequency (A/) is 
greater than or less than the local oscillator (LO) frequency. 

6. The method (33) according to claim 5 wherein step a) is further characterized by a step 
of obtaining powers O that are integer multiples of four. 

7. A detector (53) for use in a receiver (1 7) that estimates a difference (A/) between a 
received quadrature demodulated communication signal (r(t)) carrier frequency (Q and a local 
oscillator (LO) frequency used for demodulation, the detector (53) characterized by: 

a sampler (f s ) sampling the received quadrature demodulated communication signal (r(f)) 
signal components (y£t)), (y s (t)) into discrete time (y c [n])) 9 (y s [n]) signal components, said 
discrete time (y c [n]) 9 (y s [n]) signal components coupled to a complex power processor (57); 



-12- 



WO 01/82546 PCTYUS00/1 1 125 

said complex power processor (57) performs a complex power calculation on said 
discrete time (y c [n]), (y s [ri\) signal components producing real (Re{(y[n]y}) and imaginary 
(I*n{(y[n]Y}) power signal components which are output (q[n]) to a buffer (59); 

said buffer (59) accumulates a plurality (N) of complex power processor (57) outputs 
(q[n]) for a predefined period of time as a block of data (T w \ said block of data (T w ) is output to 
a complex Fourier processor (61); 

said complex Fourier transform (6 1 ) performs a complex Fourier transform on said block 
of data (T w ) outputting a frequency transform containing a plurality of frequencies (-fJ2 to (+fJ2 - 
fJN)) observed during said predefined period of time, said output (~fJ2 to (+fJ2 - fJN)) is 
coupled to a selector; and 

said selector (63) selects from said plurality of frequencies a maximum frequency which 
corresponds to the difference (A/) between the carrier (f c ) and local oscillator (LO) frequencies. 

8. The detector (53) according to claim 7, wherein said sampler (£) is further 
characterized by sampling the received quadrature demodulated continuous time (y c (t)\ (y s (t)) 
signal components at a predefined (2A/max < frequency. 

9. The detector (53) according to claim 8, wherein said complex Fourier transform (39) 
maximum frequency is further characterized with a sign (±) indicating whether the difference 
frequency (A/) is greater than or less than the local oscillator (LO) frequency. 

10. The detector (53) according to claim 9, wherein said complex power processor (57) 
is further characterized by obtaining a power O that is an integer multiple of four. 
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1 1 . The detector (53) according to claim 1 0, wherein said complex Fourier transform (3 9) 
is further characterized as a fast Fourier transform. 



12. A method (53) for estimating a difference (A/) between a received (17) quadrature 
demodulated communication signal (r(t)) carrier frequency (f c ) and a local oscillator (LO) 
frequency used for demodulation, the method (53) characterized by the steps of: 

a) sampling (Q received (r (/)) quadrature demodulated communication signal (r(f)) signal 
components (y£t)), (y£t)) into discrete time (y c [nj)), (y s [n]) signal components; 

b) performing a complex power calculation (57) on said discrete time (y c M)> (KtM) 
signal components; 

c) buffering (59) said complex power calculation output ((y[ri]Y) for a plurality (N) of 
complex power processor (57) outputs for a predefined period of time as a block of data (T w ); 

d) transforming (61) said block of data (T w ) into Fourier frequency domain values 
representing a plurality of frequencies (-fJ2 to (+fJ2 - fJN)); 

e) selecting (63) a maximum frequency from said plurality of frequencies {-fJ2 to (+/J2 - 
/,/*)); and 

f) outputting said selected frequency as the difference frequency (A/). 

13. The detector (53) according to claim 12, wherein step a) is further characterized by 
a step of sampling the received quadrature demodulated continuous time (y c (t)), (y s (f)) signal 
components at a predefined (2A/max <f s ) frequency. 
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14. The detector (53) according to claim 13, wherein step d) is further characterized by 
a step of assigning a sign (±) to said maximum frequency indicating whether the offset frequency 
(A/) is greater than or less than the local oscillator (LO) frequency. 

15. The detector (53) according to claim 14, wherein step b) is further characterized by 
a step of obtaining a power (?) that is an integer multiple of four. 

16. The detector (53) according to claim 15, wherein step d) is further characterized as 
a fast Fourier transform. 
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